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Abstract 

We measured, with the KLOE detector, the spectrum of vr+vr" invariant mass in 
a sample of 6.7 x 10^ e~^e~ 7r+7r~7 events with the photon at large polar angle 
{0^ > 45°) at a centre of mass energy y/s around the (j) mass. We observe in this 
spectrum a clear contribution from the intermediate process (p ~^ fo (980)7. ^ 
sizeable effect is also observed in the distribution of the pion forward-backward 
asymmetry. We use different theoretical models to fit the spectrum and we determine 
the fo mass and coupling constants to the (j), to tt'^tt'' and to KK. 

Key words: 13.20.-v Radiative decays of mesons, 13.20.Jf Decays of other mesons 



1 Introduction 

The 0(1020) radiative decays to fo(980) and ao(980) play an important role in 
the investigation of the controversial structure of the lighter scalar mesons [1,2]. 
At KLOE, we detect the fo through its decay to nir via the decay chain 
e~^e~ ^ — s> fo — rni'j. KLOE has already published studies on decays 
to fo7 and ao7 looking for the final states [3] and rjn^'j [4], respectively. 

On the contrary, the decay chain e+e~ ^ ^ fo(ao) KK'-f is kinemati- 
cally suppressed, and has not been observed yet. In this paper we present a 
study of the fo decay to tt+tt", based on an integrated luminosity of 350 pb~^ 
collected at the collider DA$NE during the years 2001 and 2002, at a centre 
of mass energy ^/s around the mass = 1019.45 MeV within ±0.5 MeV 
("on-peak" data). The only previous search for this decay has been published 
by the CMD-2 collaboration [5], mainly based on an energy scan around the 
mass peak. 

We look for fo n^n^ decays in events e"'"e~ n^n^^. Only a small fraction 
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of the 7r"'"7r~7 events originates from the radiative decay — *fo7 with fo 
TT+TT". The main contribution is given by e~^e~ — 7r"'"7r~7 events with a photon 
from initial state (ISR) or final state (FSR) radiation. The amplitude of each 
contribution is characterised by a different spectrum of the vr"*"??^ invariant 
mass m, and of the photon polar angle 6^ measured with respect to the beam 
axis. In particular, the ISR is the dominant contribution for small photon polar 
angles, allowing to extract the e~^e~ tt+tt" cross-section below the (j) mass 
with the so called radiative return method [6]; at large values of 6^ the ISR 
contribution is strongly reduced, so that the other processes can be observed in 
this region only. A smaller contribution comes from the decay (p p^n^ with 
— ^ 7r^7 (we call it pvr term in the following). It contributes in the low mass 
region, 400 < m < 600 MeV, with total branching ratio BR(0 p^7r^)x 
BR(p^ — ^ -7r^7)~ 4 X 10~^. Finally, the possibility to observe the decay chain 
—>fo (600)7 ~^ 7r+7r~7 is considered. 

The vr+vr" pair has different quantum numbers whether it is produced through 
FSR and fo decay or ISR: J^'-^=0"'"^ in the former case, and J^'-^=l in the 
latter. A sizeable interference term between FSR and fo decay is expected in 
the m spectrum. On the other hand any interference term between two am- 
plitudes of opposite charge conjugation gives rise to C-odd terms that change 
sign by the interchange of the two pions. Therefore, the interference between 
ISR and FSR or fo decay, results in a null contribution in the m spectrum for 
symmetric cuts on 6^ and 6'^± , and in a sizable forward-backward asymmetry, 
Ac, defined as: 

_ N{e^+ > 90°) - N{e^+ < 90°) 
' N{e^+ > 90°) + N{e^+ < 90°) ' ^ ' 

where the angle 6.,^+ is defined with respect to the direction of the incoming 
electron beam. 



2 Experimental set-up 

DA<I>NE is an e"'"e~-collider with a peak luminosity of about 10^^cm~^s~^ at a 
centre of mass energy ^/s = = 1.02 GeV. The beams collide with a crossing 
angle of 7r-0.025 rad. The KLOE detector consists of a large- volume cylindri- 
cal drift chamber [7] (3.3 m length and 2 m radius), operated with a 90% 
helium-10% isobutane gas mixture, surrounded by a sampling calorimeter [8] 
made of lead and scintillating fibres providing a solid angle coverage of 98%. 
The tracking chamber and the calorimeter are surrounded by a superconduct- 
ing coil that produces a solenoidal field 5=0.52 T. The drift chamber has a 
momentum resolution of (T{p±)/p± ~ 0.4%. Photon energies and arrival times 
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are measured by the calorimeter with resolutions oi geIE = / ^ E{GeY) 

and at = 54ps/ ^E{GeV) © 50 ps. The trigger [9] is based on the detection of 
at least two energy deposits in the calorimeter above a threshold that ranges 
between 50 and 150 MeV. The trigger includes a cosmic ray veto based on 
large energy deposits in the outermost calorimeter layers. 



3 Event selection 

We select 7r+7r^7 events by requiring a reconstructed vertex close to the in- 
teraction region with two tracks of opposite charge, emitted with polar angles 
above 45° {6j^± > 45°). We suppress the ISR component by requiring the po- 
lar angle of the total missing momentum to be larger than 45° {6^ > 45°). 
Both tracks are extrapolated to the calorimeter. A likelihood variable, based 
on the time of flight and on the shower profile (see ref. [6]), is used to select 
pions. A cut on this variable reduces the background due to e"'"e~7 events to 
a negligible level. 

In order to remove 7r'^7r~7r° and /i"'"/i~7 events, we define the track mass vari- 
able Mt as the solution of the equation: 

|P0 -Pi -P2I =E^- ^pj + Ml - + Ml, (2) 

where pi and p2 are the momenta of the two tracks, and where E^ and p^ are 
the (j) energy and momentum, respectively. These are evaluated run by run 
using samples of Bhabha scattering events. Eq. (2) is verified by events with 
two particles of mass Mt and a third particle with null mass. Mt is required 
to be equal to the pion mass within ±10 MeV . To reduce the residual 7r"'"7r~7r° 
contamination and to remove badly reconstructed e~^e~ — > vr'^vr" events with 
soft photon emitted, we require a calorimeter cluster matching the missing 
energy and momentum. The cluster is required to be non associated to tracks, 
to have an energy E^ > 10 MeV, and to have a time compatible with a pho- 
ton coming from the interaction vertex, — R^/c\ < 5at{E^), where is the 
cluster time, is the flight distance, and at{E^) is the time resolution for 
photons of energy E^. The requirement E^ > 10 MeV translates in an effec- 
tive cut m < 1009 MeV. Finally, the angle Q between the missing momentum 
and the photon direction derived from the cluster position, has to be below 
0.03+3/-E^(MeV) rad. The dependence of the Q cut on E^ reflects that of the 
cluster position resolution on the photon energy. 

We select 6.7x10^ events. The spectrum of the vrvr invariant mass m for these 
events and the forward-backward asymmetry dependence on m are shown in 
Fig. 1. We observe a small bump in the m spectrum in the region where the 
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fo(980) is expected to be. A signal is observed also as a dip in the forward- 
backward asymmetry for the same values of m. 




500 600 700 



900 1000 
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900 1000 
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Fig. 1. (a) vr'^vr" invariant mass spectrum of the selected sample. The spectrum is 
dominated by the ISR component, showing the p-co interference pattern. The signal 
of the fo (980) appears as a small peak around 980 MeV. The drop for m > 1000 MeV 
is due to the drop of the detection efficiency for low energy photons, (h) For- 
ward-backward asymmetry defined in Eq. (1) as a function of m. The dip in the 
region of the fo (980) is evident. 



Total efficiency and residual background distributions are shown in Fig. 2, 
as evaluated by Monte Carlo with corrections based on data control samples 
[10]. The simulation of the ISR and FSR contributions is based on the EVA 
generator [11]. The efficiency decrease at low masses is due to the increased 
occurrence of low-p^ pions with 9.„± < 45° that escape the selection; the 
decrease for higher masses, starting from ~ 800 MeV, is partly due to the 
cosmic ray veto and partly to the photon detection efficiency. In fact, high 
momentum pions, which deposit large energy in the outermost calorimeter 
layers, veto the event with high probability. Moreover, low energy photons 
{E.y < 20 MeV) are detected with an efficiency lower than 80%. The efficiency 
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Fig. 2. (a) Total efficiency as a function of m. (b) The Monte Carlo expected 
contributions of the main background sources, vr+vr^vr'^ (open circles) and fj,'^fj,~^ 
(crosses), normalised to the integrated luminosity, and compared to the data spec- 
trum. 



for the cosmic veto is evaluated using samples of pre-scaled events with no 
veto applied. The photon detection efficiency is measured as a function of 
from 7r"'"7r~7r° and e"^e~7 control samples [8]. 

After the selection, — n^n^ir^ decays give the only significant contribution 
to the background. 



4 Description of the fit 



We fit the TT+TT invariant mass spectrum, dN/dm, with the function: 



— = L- te{m) ( ^^^^^ + ^^^^^ + ^^P^ _^ ^^scal _^ C?q-4al+FSR \ _^ 

dm \ dm dm dm dm dm J 

back, (3) 
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where L^t is the integrated luminosity, e(m) is the selection efficiency, and 
back is the residual background. The ffist three terms in parenthesis are here 
called the "non-scalar" terms. The analytic expressions for the first and sec- 
ond terms, ISR and FSR, are taken from ref. [12], while the pvr term is taken 
from ref. [13]. The pion form factor [14], entering the ISR term, depends on 
the masses and widths of the p", u and p' mesons, and on the two non di- 
mensional parameters a and (3, which correspond to the sizes of the u and p' 
contributions, respectively. We leave the quantities mpo, TpO, a, and (3 as free 
parameters of the fit while the masses and the widths of the u and p' mesons 
are fixed to the PDG values [15]. The pvr term is multiplied by a scale factor, 
a^TT, which is expected to be equal to unity. If ap,r=l the number of pn events 
corresponds to approximately 1% of the total, broadly distributed in the low 
mass region. The possible interference between the pn and the scalar terms 
is neglected. The last two terms in parenthesis, seal and scal+FSR, depend 
on the amplitude for the decay — >fo7 — ^ 7r"^7r~7, the latter being the inter- 
ference term between fo and FSR. The scal+FSR term can be either added 
(constructive interference) or subtracted (destructive interference). 



We perform three fits corresponding to three different approaches in the de- 
scription of the scalar amplitude. 

The first fit is the Kaon-Loop fit (KL) [1,12]: the couples to the scalar 
through a loop of charged kaons. The formalism allows the inclusion of more 
than one scalar meson. For each scalar meson there are three free parameters 
of the fit: the mass and the couplings to K^K~ and to tt+tt^. For the fo scalar 
meson only, the amplitude reduces to: 

^x.=^(m^)e-^^-) f°^"^;^;;:;-- , (4) 



where s is the square of the centre of mass energy, g{gK+K+ and gign+w- are the 
two couplings, g{m?) is the kaon-loop function [1], 5{m) is the phase of the 
background to the tttt elastic scattering, and is the fo inverse propagator 
with the finite width corrections [1]. 

In the second fit, called No-Structure fit (NS) [16], a direct coupling 5'(/,fo7 of 
the to the fo is assumed, with a subsequent coupling gi^^T,+T,- of the fo to 
the 7r"'^7r" pair. The fo amplitude is a Breit-Wigner with a mass dependent 
width [17] added to a polynomial complex function, the continuum, to allow 
an appropriate dumping of the resulting line shape. The amplitude depends 
on eight parameters: the mass mf^,; the three couplings g^i^-y, gioir+n-, and 
gfoK+K-] four parameters describing the continuum: two coefficients qq and ai 
and two phases and bi. The amplitude is: 

^ ^#o7^fo.+^- ^ ^^.6oP.M + ai "''"'"'^ '^^^-^-), (5) 
Ao("^) mi ml 
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where p-Kim) = y 771^/4 — is the vr momentum in the fo rest frame and 
Di^{m) is the fo inverse propagator: 



where g^,,^^ = ^J^gi,,T,+-K- and q^^kk = 9ioK+K- = 9{^koko- These couphngs 
have the same meaning as those defined within the KL frame and are related 
to the non dimensional couplings (7^ and (see for instance [18,19]) through 
the relations g^ ~ dlnn/ i^^'f^D and gx ~ sIkk/ ('^'^'"^l) strictly valid only 
when m ~ mfp. In order to obtain the correct phase behaviour consistently 
with chiral perturbation theory predictions [20], bo is expressed as a function 
of the other parameters, reducing the free parameters to seven. 



Ami 



+ 



9loKK 

IGvr 



Arn\± 



(6) 



Finally in the Scattering Amplitudes (SA) fit [21] the amplitude is the sum of 
the scattering amplitudes Tn = T^-k-k nir) and T12 = T(7r7r KK), whose 
shapes are fixed by independent experimental information [22]: 



Asa = {m- 171^) 

(02 + hrn^ + C2ni^)Ti2 




(ai + him^ + Cim^)Tii + 

A 



(7) 



free parameters are the six coefficients of the two polynomials, mg and the 
phase A. Once the amplitude is determined by the fit, it is analytically con- 
tinued in the complex m plane and the coupling is determined by the pole 
residue. The coupling gf^, having the dimension of an energy, is connected to 
the partial width r(</) — > 7fo(980)) through the relation [21]: 



7fo(980)) = 



5 Results 



We fit the data in the region 420 < m < 1010 MeV, using bins 1.2 MeV wide 
[23]. 

First we discuss the fits KL and NS. In both destructive interference 
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Fig. 3. Result of the KL fit (a)-(b)-(c) and of the NS fit (d)-(e)-(f). (a)-(d) Data 
spectrum compared with the fitting function ( upper curve following the data points ) 
and with the estimated non-scalar part of the function (lower curve); (h)-(e) fit 
residuals as a function ofm; (c)-(f) the fitting function is compared to the spectrum 
obtained subtracting to the measured data the non-scalar part of the function in the 
fo region. 

is preferred by the fit for the {da / d'm)l^^_^p^j^ term, while the constructive 
interference is strongly disfavoured. The results are shown in Fig. 3, the of 
the fits and the values of the parameters are given in Tab.l. The non-scalar 
part is well described by the parametrisation used, while we are clearly not 
sensitive to the pn term. The fo signal appears as an excess of events in the 
region between 900 and 1000 MeV. In the KL fit, the attempt to include a 
second scalar meson (the fo(600) or a), with either E791 [24] or BES [25] 
masses (respectively 478 and 541 MeV), and with free couplings, gives no 
improvement to the fit. Moreover, since the couplings preferred by the fit are 
compatible with zero, within the statistical errors, the fo(600) is unnecessary 
to describe these data. 

After the subtraction of the non-scalar part obtained in the KL and NS fits. 
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an asymmetric peak around 980 MeV with a FWHM of 30— 35 MeV and a 
height ~25% of the total is obtained, as shown in Fig. 3(c) and (f). Such a 
peak does not directly represent the fo shape but it results from the sum of a 
broad term {da/dm)scai and a negative interference term ((i(T/(im)scai+FSR that 
cancels the low mass tail. The NS fit requires a significantly larger value of /5 
than the KL fit does. This results in a non-scalar part ~ 4% larger in the fo 
peak region, and hence a correspondingly smaller signal size. 

Table 1 

Parameter results andx^ of the two fits KL (kaon-loop) and NS (no- structure). The 
results given in parentheses are not directly parameters of the fits hut are evaluated 
as functions of the fit parameters. 





KL 


NS 


(p(x')) 


538/483 (4.2%) 


533/479 (4.4%) 


mf(, (MeV) 


983.0±0.6 


977.3±0.9 


^#07 (GeV-^) 




1.48±0.06 


9ioK+K- (GeV) 


5.89±0.14 


1.73±0.12 


9fon+n- (GeV) 


(3.6) 


0.99±0.02 




2.66±0.10 


(3.1) 


ao 




6.00±0.02 


ai 




4.10±0.04 


6i(rad/GeV) 




3.13±0.05 


mpo (MeV) 


773.1±0.2 


773.0±0.1 


r^o (MeV) 


144.0±0.3 


145.1±0.1 


a (xl0~3) 


1.65±0.05 


1.64±0.04 


/3 (xl0~3) 


-123±1 


-137±1 




O.OibO.6 


1.5±1.4 



In order to assign the uncertainty to the parameters extracted from the fits, 
we have done a study of several systematic effects [23] . We have repeated the 
fits by varying the following quantities: the luminosity value around its to- 
tal estimated error (~ 1% [26]); the shape of the photon efficiency and the 
linearity of the photon response curves; the size of the ir'^iT'ir^ background; 
the bin size, and the start and end points of the fit. While repeating the fits, 
the parameters of the non-scalar part are held fixed to their baseline values. 
Finally in order to take into account the systematic effect due to the limited 
knowledge of the non-scalar part of the spectrum, the NS fit has been repeated 
using the non-scalar parameters obtained from the KL fit and vice-versa. In 
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Table 2 

Intervals of maximal variations for the /o parameters resulting from the systematic 
uncertainties studies done on both fits. Notice that the intervals obtained are larger 
than the fit uncertainties given in the previous table. 



parameter KL NS 



(MeV) 


980- 


987 


973- 


981 


5foK+K- (GeV) 


5.0 - 


6.3 


1.6- 


2.3 


ftoTT+TT- (GeV) 


3.0- 


4.2 


0.9- 


1.1 




2.2- 


2.8 


2.6- 


4.4 


5#o7 (GeV-^) 






1.2- 


2.0 



Table 3 

Results of SA fit. The o,nd the numerical values of the parameters are given 
together with the resulting values for the parameters of the non-scalar part. 



(p(x')) 


577/477 (0.1%) 








ai 


11.9 


a2 


-14.7 


nipo (MeV) 


774.4±0.2 


bi 


3.3 


b2 


-15.3 


r^o (MeV) 


142.8±0.3 


Cl 


-15.1 


C2 


35.8 


a (xlO^^) 


1.74±0.05 


mo 


0. 


A 


-1.63 


P (xlO-3) 


-100±18 












0±2 



Tab. 2 we give the maximal variation intervals for the parameters, resulting 
from the studies discussed above. 

The two fits have slightly overlapping intervals for the fo mass, and are both 
in agreement with the PDG interval 980± 10 MeV. We observe a large dis- 
crepancy between the KL and NS couphngs gioTr+Tv- and gt^^K+K-- The KL 
fit gives couplings in reasonable agreement with the KLOE results obtained 
with the final state 7r'^7r''7 [3]. The two fits are in agreement on the ratio 
R = 5ff^j^+j^_/(yf|j^+^_, pointing to an fo more coupled to kaons than to pi- 
ous. Finally if we define an effective branching ratio as the integral over 
the full spectrum of the fo term normalised to the total width, we obtain 
BR(0 ^fo(980)7)x BR(fo(980)^ n+n^)= 2.1x10"^ and 2.4x10^^ for KL 
and NS fits respectively. 

The SA fit is shown in Fig. 4 and in Tab. 3. The is poorer especially in 
the fo peak region. Notice that a much better can be hardly expected since 
Til and Ti2 are derived from data sets that are less accurate than the data 
presented here. In any case by properly normalising the amplitude we obtain a 
value ~ 6.6 x 10~^GeV, in agreement with the value obtained in ref. [21] by 
fitting the KLOE data on TT^ir^'y together with other data. However we stress 
that such a value corresponds to an effective branching ratio of ~ 3 x 10"^, one 
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Fig. 4. Same as Fig. 3 for fit SA. In this case the fit lowers the non-scalar part and 
gives a larger signal than in the KL and NS fits. In the peak region the fit is clearly 
poorer than KL and NS fits. 

order of magnitude lower than the one obtained from the other two fits. This 
can be understood since in this case the fit prefers a constructive interference 
term, hence the scalar term has a smaller size. 



Using the results of the KL fit, we predict the dependence of the cross-section 
a(e+e- ^ 7r+7r-7,45° < 6^ < 135°, 900 < m < 1000 MeV) on yfs. The pre- 
dicted behaviour is compared to the data in Fig. 5, in the y/s range between 
1016 and 1023 MeV. Besides the "on-peak" data sliced in 0.1 MeV wide bins, 
we show two "off-peak" points taken at ^/s = 1017 and 1022 MeV, respec- 
tively. We observe a good agreement for the on-peak data, and a marginal 
agreement for the two off-peak points. 
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Fig. 5. Centre of mass energy dependence of the cross-section for events with m 
in the range 900 - 1000 MeV. The open points are the "on-peak" data sliced in 
0.1 MeV wide bins, the full points are the ''off-peak" data. The curve is the absolute 
prediction based on the KL fit parameters. 




Fig. 6. The forward-backward asymmetry data (full circles) compared to the Monte 
Carlo expectations based on the non-scalar part of the spectrum only (open tri- 
angles), and on the non- scalar plus /o part obtained from the KL amplitude (open 
squares). The right plot shows the detail of the comparison in the fo region. 

Finally, following the suggestion contained in ref. [27], we compared the be- 
haviour of the forward-backward asymmetry as a function of m, shown in 
Fig. 1(b), with a simulation including the fo contribution besides the ISR and 
FSR [28], and the effect of the 7r"'"7r~7r° background that dilutes the asym- 
metry in the low mass region. The comparison is shown in Fig. 6. The KL 
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parametrisation of the fo amplitude has been used. The inclusion of the fo7 
term is essential to have an acceptable agreement between data and simulation 
in the region of the fo peak and also in the low mass region. This means that 
the low mass tail of the fo7 amplitude that is cancelled in the m spectrum by 
the destructive interference with FSR, is on the contrary well evident in the 
Ac spectrum due to the interference with ISR. 



6 Conclusion 

Summarising, we found a clear evidence for the process —^fo (980)7 ~^ 7r+7r~7 
in the vrvr invariant mass spectrum and in the behaviour of the forward- 
backward asymmetry. An acceptable description of the data is obtained with 
fits KL and NS. Both fits predict the fo to be strongly coupled to kaons, the 
ratio R between g'j^x+K- ^^"^ fi'foTr+vr- being well above 2 as in the KLOE mea- 
surement of fo 7r°7r° [3]. The coupling to the 0, f^^fo^ is found using the 
NS approach and is in the range 1.2 — 2 GeV~^. A marginal agreement is 
obtained by applying the fit SA. 
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